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Abstract 

Energetic particles in a turbulent medium can be subject to second- 
order Fermi acceleration due to scattering on moving plasma waves. 
This mechanism leads to growing particle momentum dispersion and, 
at the same time, increases the mean particle energy. In the most fre- 
quently met situations both processes can be represented by a single mo- 
mentum diffusion term in the particle kinetic equation. In the present 
paper we discuss the conditions allowing the additional term for regular 
acceleration to arise. For forward-backward asymmetric scattering cen- 
tres, besides the diffusive term one should explicitly consider the regular 
acceleration term in momentum space, which can consist of the first- 
order (oc V), as well as the second-order (oc V 2 ) part in the wave velocity 
V. We derive the condition for the scattering probability in the wave 
rest frame required for vanishing the regular acceleration term and pro- 
vide a simple mechanical example illustrating the theoretical concepts. 
Finally, we address its possible role in cosmic ray acceleration processes. 
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1 Introduction 

The original theory conceived by Fermi [5] for the acceleration of cosmic rays 
considers charged particles scattered by magnetic clouds moving with speeds 
V relative to the plasma rest frame. The resultant acceleration is represented 
by a systematic gain term. The following studies considered some particular 
acceleration scenarios involving large amplitude magnetohydrodynamic waves 
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and allowing for regular acceleration of cosmic ray particles (cf. Parker [7]). A 
general approach is based on the Fokker-Planck equation (cf. Parker & Tidman 
[8]; Clemmow & Dougherty [2]; Lacombe [6]; Blandford & Eichler [1]). In the 
recent years this approach supplemented with the quasi-linear derivation of 
transport coefficients for cosmic rays interacting with Alfvenic turbulence (e.g. 
Schlickeiser [9]) became an analytical basis for considering energetic particle 
transport phenomena in rarefied astrophysical plasmas. Within this approach 
the detailed balance principle applied to particle interaction with the field of 
magnetized plasma waves leads to the reduction of the general Fokker-Planck 
equation to the pure momentum diffusion equation. Let us note, however, that 
the principle of detailed balance which arises in a natural way in the kinetic 
description of gas of scattering particles is no longer a natural assumption for 
particles scattered on external 'heavy' scattering centres. In our opinion this 
till now neglected fact deserves more detailed study because of its potential 
importance for cosmic ray particle acceleration in turbulent media. 

In the present work we consider detailed conditions for the angular scatter- 
ing probability that is required for reduction of the Fokker-Planck equation de- 
scribing particle momentum scattering into the momentum diffusion equation. 
The next section presents general derivation leading to such simplification. 
Then, in section 3, we present calculation of the Fokker-Planck coefficients for 
the scattering centres (below we call them simply 'waves') moving with ve- 
locities much smaller than the particle velocity. We show that the scattering 
probability distribution, which is - in the scatterig wave rest frame - symmetric 
with respect to the interchange of particle momentum before and after scat- 
tering, is required to describe the process as a pure momentum diffusion. A 
simple mechanical model allowing the breaking this symmetry and yielding a 
regular acceleration (deceleration) term is discussed in section 4. In the dis- 
cussion (section 5) we stress the importance of the possibility to have such 
regular term in the stochastic acceleration process for cosmic ray production 
in astrophysical sources. 

2 The Fokker-Planck approach 

As above mentioned the approach based on the Fokker-Planck equation is 
generally used to describe energetic particle interactions with magnetohydro- 
dynamic turbulence. Following Blandford & Eichler [1] we write — in the 
phase space — the equation describing the evolution of the isotropic part of 
particle distribution function /(r, p, t), due to momentum scattering as 
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dt Vp At J + 2 Vp 



(ApAp) . 
At J 



(2.1) 



where the Fokker-Planck coefficients, the mean momentum gain rate rep- 
resents the Fermi term and the rate for growing momentum dispersion is gov- 
erned by the term ^ A ^ p ^ ; V p = dj 'dp. With a probability for changing a 
particle momentum in an individual scattering from p to p' (= p + Ap) given 
by the function ^(p',p) , with the normalization / ^dAp = 1, these coeffi- 
cients are given as 



(A P ) 
At 

and 



(ApAp) 



^/*(p',p)(p'-p)dAp (2.2) 

At ^/*( P ',p)(p'-p) 2 rfAp , (2.3) 

where 1/At is the scattering frequency. If the principle of detailed balance 
holds, tf(p',p) = *(p,p'), we have (Blandford & Eichler [1]): 

The equilibrium momentum distribution can exist on the condition of the 
momentum-change velocity vanishing (i.e. a constant vector in parentheses in 
Eq. 2.4). Then, the Fokker-Planck coefficients relate as 

(Ap) 1 (ApAp) 

At "2 Vp At ' (2 ' 5j 

and Eq. (2.1) takes a simple form 

^ = V p • {D pp ■ V p /} , (2.6) 

where one introduced the momentum diffusion tensor D pp = ( Ap Ap) / (2 At) . 
For the isotropic momentum diffusion D pp = D p I, where D p is the momen- 
tum diffusion coefficient and I is the diagonal unit matrix. Then Eq. (2.1) 
transforms into the form with the scalar momentum and the operator dj 'dp 
replacing the vector p and V p , respectively, and, accordingly, Eq. (2.6) takes 
the form 
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df_ _ }_d_ \ 2 D 9f\ 

dt p 2 dp\^ p dpj 

In the present work we consider detailed conditions required for reduction of 
Eq. (2.1) into (2.7). 



3 Fokker-Planck coefficients and the condition 
for vanishing a regular acceleration term 

Any wave damping or boosting process in rarefied astrophysical plasmas can 
be usually neglected for MHD waves at time scales relevant for single act of 
energetic particle interaction with a given wave. Due to this fact the elec- 
tric field can be assumed to vanish in the wave rest frame ('wave frame') and 
the transport coefficients in the momentum space can be obtained using a 
straightforward procedure. Let us consider an individual wave (scattering cen- 
tre) taken from the ensemble of waves with isotropic distribution of velocities. 
In the wave frame the scattering conserves the particle energy and can be 
described with the angular scattering cross-section. Thus, the energy gain of 
any individual particle can be derived from two Lorentz transformations, the 
transformation from the mean plasma rest frame ('plasma frame') to the wave 
frame before the scattering, and in the opposite direction after the scatter- 
ing. To derive transport coefficients one simply has to calculate the change in 
momentum on scattering and to average it over incoming and outgoing pitch 
angles. This approach was proposed by Duffy [4] and the discussion in this 
section follows his work in essentials. 

Let us take the z-axis of the considered reference frame along the mean 
background magnetic field. In the plasma frame we denote an X quantity 
before interaction and after interaction with a prime, X' . In the wave frame 
an index w is added to respective quantities. Then, we denote with V a 
wave velocity assumed to be along the magnetic fieldf], with v being a particle 
velocity (v ^> V), and with p and fi a particle momentum and a pitch angle 
cosine, respectively (below, \x is called shortly 'a pitch angle'). The Lorentz 
transformation of particle momentum into the wave frame gives, to second 
order in V/v (as required by the Fokker-Planck approach), 

1 But the derivations are performed with all angles measured with respect to the direction 
of the scattering centre velocity and, in principle, has nothing to do with the magnetic field. 
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Pw=P 



fiV 1 , A V 2 fi 2 V 2 
v 2 v 1 v z 2 c z 



and, to the first order in V/v, 



V 



fi w = fi - (l - fi 2 ) - 



(3.1) 



(3.2) 



The particle momentum change Ap = p' — p at an individual scattering trans- 
forming \i into fi' is 



Ap = 



/2 



/i 



2^ 



2 c 2 



P 



(3.3) 



where the terms of the higher-than-second order in V/v are neglected. Let 
us note that in the expression (3.3), there are no zero-order terms in wave 
velocity. Therefore, in order to derive < Ap > and < {Ap) 2 > valid up to 
the second-order in V/v one should derive the required weighting factors up 
to the first-order and to the zero-order, respectively. The factor which bears 
all physical information about the scattering is a probability for the particle 
pitch angle fi to be transformed into fi'. The scattering process is described 
here in the wave frame, where we consider the probability P w (p w , fi w , fi' w ) for 
a particle with momentum p w to change fi w to fi' w during a time At w (= At + 
terms 0(V 2 /v 2 ) ). At a scattering the particle momentum p w is assumed to 
be conserved. To the first order the scattering frequency in the plasma frame 
is proportional to the Lorentz factor 



7 v 



V 



(3.6) 



where 7 = — v 2 . Thus, expressions for the Fokker-Planck coefficients 

valid to the second order in V/v are 



(Ap) = 1 £ dfi £ dfi' - fj^j P(p, fi, fi') Ap(p, fi, fi') (3.8) 

and 

((Ap) 2 ) = i £ dfi £ dfi' P(p, fi, fi') Ap 2 (p, fi, fi') , (3.9) 
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where the probability P should be expressed with the use of P w , by expanding 
it to the first-order in V/v. Let us note that the above coefficients are the ones 
for the isotropic part of the distribution function, denoted above as f(p). In 
the wave frame, a number of particle transitions (/i w — > ji' w ) per unit time is 
proportional to P w (p w , /i w , fx' w ) and to the slightly anisotropic particle density 
fw(Pw, fJ> w ) — /(p) - (d///d/x w ). Therefore, one can write the probability function 
in the plasma frame as 

P(p, n, //) = P w [p w {p, n), fi w {fi),fi' w {fi')} -^f . (3.10) 

With the use of Eq-s (3.1,2) the integrand in (3.8) can be expanded to the 
first-order as 



1 v \ m >\ P f v ,v\ v dpj^y) 

1 - //— P(p, fi, ii ) = P w (p, fi, fi ) 1 - fi— + 2/i — - pfj, 

V J \ V v J v op 



(l n 2 ) - 9Pw (l n' 2 ) - 9Pw ^ ^ ^ (3 11) 

In the integral in (3.9) we need analogous expansion to the zeroth-order in 
V/v, P{p,n,n') = P w (p, /i, //), to obtain the momentum diffusion coefficient 

as 



Dp = ^ Q 2 /-i d/i W - A* P»(P, V, *t) ■ (3-12) 

In the above formula one can consider the transition probability in the wave 
rest frame, where the cosmic ray scattering can be described in the most simple 
way. For example, in the case of scattering on magnetohydrodynamic waves 
the electric field vanishes in such a frame leading to interactions conserving 
particle energies. 

Let us consider the condition (2.5) required to reduce the acceleration pro- 
cess to the purely diffusive one (Eq. 2.6,8). With the scattering probability 
given in Eq. (3.11), one can derive the mean momentum change in Eq. (3.8). 
For a general form of P w the regular acceleration term should be appended to 
Eq. (2.7) besides the momentum diffusion term. Then, it takes a more general 
form 
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dt p 2 dp pp dp J p 2 dp ^ ^ Teg > 



(3.13) 



where p reg = < Ap > reg /At is a regular momentum changing term, equivalent 
to the regular acceleration if p reg > 0. With Eq-s (2.5,3.8-9) one obtains 



< Ap > 



reg 



1 fV\ r 1 , r 1 



p 4 V v 

1 fV\ 2 r 1 . r 1 



J ^ dfj, J fy' (a*' - A-0 ^to(p, A 1 ? A*') 



~4 \~v) J-i d ^ i-i _ ^) /i ' ^ 



+ 8 p Uj 7-1 ^ 7-1 ^(^-^ 



dP w (p, fi, n') 
dp 



(3.14) 



The analytic condition for vanishing of the term (3.14) is the symmetry of the 
scattering probability P w with respect to the interchange of \i and // 



P w (p, fi, //) = P w (p, fi' , fi) . (3.15) 

If the principle of detailed balance is satisfied, the above condition must also 
held (in general, the opposite is not always true). Then, for two opposite 
scattering acts (p — > p') and (p' — > p) the respective momentum gain and 
loss terms cancel each other exactly and any net energy change comes from 
different particle densities near p and p'. In most cases of ordinary scattering 
processes the condition (3.15) is trivially satisfied and the momentum diffusion 
fully describes the acceleration process (cf. Drury [3]). In order to illustrate the 
physical reasons for non- vanishing regular term p reg , in the following section a 
simple mechanical example with asymmetric particle scattering is presented. 



4 An example with a non-vanishing regular 
acceleration term 

As a simple illustration for the scattering process with the cross-section show- 
ing the forward-backward asymmetry we consider the elastically reflecting reg- 
ular tetrahedrons (i.e. triangular pyramids with all 4 walls being the equal 
triangles), moving always with one of its bases directed forward and perpen- 
dicular to its direction of motion. Let us measure all angles with respect to the 
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base normal, pointing by definition into p — 1. We consider particle reflections 
in the tetrahedron rest frame. If we take a stream of particles hitting the base 
with p — — 1, then the reflected particles will be characterized with pitch angle 
p' — 1. However, if a stream hitting the tetrahedron will have originally p = 1, 
then no one of particles reflected from the inclined walls of tetrahedron will 
be scattered into pi = —1. This simple example demonstrates the breaking of 
the condition (3.15). Next, let us derive the diffusion coefficient and, both, the 
diffusive and the regular acceleration term existing for such scattering centres. 
From numerical integration of Eq-s (3.12) and (3.14), with the reflection prob- 
ability defined by the tetrahedron geometric shape, we obtain the following 
transport coefficients 



where we put At — 1. In expressions (3.14) and (4.3) for p reg one may note 
the existence of terms representing both the first-order and the second-order 
regular acceleration. The occurrence of the term linear in V/v is of particu- 
lar interest in any application of the theory to particle acceleration. In our 
example, for small V/v, this term can dominate over another terms. 

5 Final remarks 

In the present work we considered conditions required for reduction of the 
stochastic particle momentum scattering process into a pure momentum diffu- 
sion. Our considerations were intended to translate this rather general state- 
ment into the terms of scattering probability in particle pitch angle and trans- 
port coefficients derived to the second-order in V/v, as required by the Fokker- 
Planck approach. The derived formula requires an equal probability in the 
wave frame to scatter a particle from the pitch angle p to pi and from pi to p. 
It is equivalent, after averaging over the particle phase angle, to assumption 
of the detailed balance principle. We also derived an expression for a regular 
acceleration term occurring when the symmetry of the above probability is bro- 
ken. In the previous section we illustrate with a simple example the conditions 
allowing for the regular acceleration term to exist. The resulting coefficient 




(4.1) 




(4.2) 




(4.3) 
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p reg can consist of the linear in the scattering centre velocity term as well as 
the second-order term proportional to (V/v ) 2 . In astrophysical conditions one 
usually has 1/ < c. Thus the non-vanishing linear term could provide a fast 
way to accelerate particles in the first-order process, with the regular term up 
to three orders of magnitude larger than the diffusive term. 

The importance of the present derivation depends on the possibility to oc- 
cur conditions allowing to break the scattering symmetry in the wave frame 
(Ostrowski & Siemieniec-Ozi§blo, in preparation). In this respect one can 
consider regions of space with MHD turbulence influenced by the strong ex- 
ternal driving force and leading to high-amplitude waves. Thus, particles can 
be scattered by non-linear asymmetric waves, allowing for different particle 
distributions resulting from head-on and tail-on collisions. We expect such 
conditions to occur near high Mach number shock fronts as well as in bound- 
ary regions of high velocity jets penetrating through the cold medium. 
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